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1.  Introduction 


The  characterization  of  the  scattering  responses  of  nonlinear  circuit-element-loaded  scatterers  is 
an  important  problem  in  both  electromagnetic  interference  and  remote  sensing  applications.  As 
pertinent  to  communication  systems  design,  it  is  well-known  that  nonlinearly  modulated 
scattering  phenomena  constitute  a  factor  in  understanding  how  inter-coupling  effects  among 
different  radios  influence  transceiver  perfonnance  and  signal  fidelity.  Fundamentally,  from  the 
sensing  perspective,  a  nonlinearly  loaded  target — when  properly  excited — could  potentially  re¬ 
radiate  unique  features,  enabling  it  to  be  differentiated  from  an  ordinary  scatterer.  For  example, 
the  antenna,  filter,  and  amplifier  block  at  the  RF  front-end  of  the  receiver  of  an  electronic  device 
can  modulate  a  probing  signal,  facilitating  opportunities  for  signature  classification  and 
detection — in  particular,  the  presence  of  nonlinear  components  can  lead  to  hannonic  signal 
generation  within  the  electronic  system  and  subsequent  unintended  re-emissions  through  the 
antenna  port;  these  re-emissions,  in  turn,  can  be  exploited  to  aid  the  extraction  of  the  target 
embedded  in  linear  clutter.  From  the  radar  implementation  perspective,  the  tantalizing  prospect 
of  "clutter- free"  detection  has  been  previously  discussed  in  literature  in  conjunction  with  the 
introduction  of  the  hannonic  radar.  ’  Within  the  above  prescribed  context,  the  electromagnetic 
characterization  of  nonlinearly  loaded  antennas  in  itself  is  seen  to  be  an  important  sub-problem 
that  deserves  further  in-depth  consideration.  Accordingly,  in  this  study,  both  the  harmonic 
scattering  and  imaging  responses  of  nonlinearly  loaded  antennas  located  in  the  presence  of  a 
half-space  environment  or  within  a  complex-room  environment  are  investigated — with  specific 
relevance  to  active  RF  sensing,  detection,  and  localization. 

The  electromagnetic  behavior  of  nonlinearly  loaded  radiators  has  been  a  subject  of  interest  in 
many  studies.  A  variety  of  direct  time-domain  integral  equation  methods  have  been  proposed  for 
analyzing  nonlinearly  loaded  wire  antennas.3  6  Volterra  series  analysis  is  adapted  in7  for 
modeling  the  nonlinear  antenna  system  after  transforming  the  electromagnetic  problem  into  a 
nonlinear  network  using  the  method-of-moments.  The  focus  of  more  recent  studies  has  been  on 
frequency  domain  methods:  in,8  a  hannonic  balance  technique  is  introduced  for  simulating 
scattering  from  wire  scatterers  with  nonlinear  loads;  in,9  12  several  techniques  are  proposed  for 
treating  either  the  single  dipole  or  an  array  of  dipoles,  including  the  reflection  method,  the 
application  of  neural  networks,  and  the  harmonic  balance  method  in  combination  with  the 
inexact  Newton  approach  or  the  particle  swarm  algorithm.  All  of  the  aforementioned  studies 
assume  the  antenna  is  in  a  free  space  environment;  the  inclusion  of  lossy  ground  effects  in 
calculating  the  nonlinear  response,  however,  is  considered  in13  16 — where  the  antenna  is  placed 
above  the  ground  surface.  It  also  should  be  mentioned  that  the  scopes  of  these  studies  do  not 
fully  encompass  the  remote  sensing  problem — that  is,  the  discussions  therein  do  not  consider  the 
localization  of  the  loaded  radiator  by  exploiting  the  nonlinear  scattering  response. 
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Spatial  and  spectral  array  signal  processing  techniques  can  be  used  for  interpreting  the 
information  embedded  in  the  scattering  returns  in  order  to  localize  the  nonlinearly  loaded 
antenna  target.  For  example,  in  active  sensing — analogous  to  the  case  of  the  ultra-wideband 
linear  radar  system — after  sweeping  over  a  band  of  interrogation  frequencies  and  subsequently 
capturing  the  scattering  responses  at  the  hannonic  frequencies,  an  extension  of  the  conventional 
time-reversal  method17  19  may  be  employed  for  target  imaging  purposes.  The  application  of  such 
a  sensing  approach  would  enable  the  target  to  be  isolated  from  surrounding  linear  clutter; 
however,  a  rigorous  study  in  this  area  has  not  been  carried  out  in  existing  works.  A  direction-of- 
arrival  estimation  technique — or  its  many  variants,  as  applied  to  the  hannonic  returns — can  also 
be  used  to  resolve  the  angular  position  of  the  nonlinear  source;  furthermore,  direction-of-arrival 
and  time-of-arrival  processing  can  be  combined  to  obtain  the  absolute  position  of  the  target.  A 
survey  of  the  literature  shows  that  in  the  harmonic  radar  community,  several  approaches  are 
currently  exploited  for  target  tracking:  for  instance,  in,“  a  harmonic  radar  direction  finder  is 
discussed;  and  in,“  a  modified  frequency  modulated  continuous  wave  radar  is  proposed  for 
target  range  determination. 

In  this  work,  an  investigation  is  performed  to  explore  the  feasibility  and  the  mechanisms  of  using 
a  subspace-based  method  for  imaging  a  scene  consisting  of  nonlinearly  loaded  antenna  targets 
and  linear  clutter.  The  focus  here  is  primarily  on  the  half-space  problem,  as  the  study  is  intended 
to  facilitate  the  development  of  a  complementary  nonlinear  detection  capability  for  an  existing 
near-ground,  forward-looking  linear  imaging  radar  system  introduced  in.  As  such,  the 
theoretical  formalism  is  first  established  purely  in  terms  of  the  half-space  detection  problem 
before  the  complex-room  sensing  problem  is  considered  as  a  special  case.  The  overall  work  is 
put  forth  as  follows:  in  Section  2,  the  analysis  and  simulation  framework  is  delineated — first  the 
formalisms  of  the  solver  for  analyzing  the  linear  and  nonlinear  multistatic  scattered  fields  from 
nonlinearly  loaded  antennas  are  outlined,  then  the  subspace  imaging  algorithm  for  processing  the 
hannonic  re-radiation  is  presented;  in  Section  3,  the  electromagnetic  scattering  and  propagation 
solver  and  the  harmonic  imaging  technique  are  applied  to  study  the  detection  and  localization  of 
nonlinear  targets  in  half-space  and  complex-room  environments;  and  in  Section  4,  a  summary  of 
the  investigation  is  given. 


2.  Analysis  and  Simulation  Framework 


As  a  starting  point,  standoff  sensing  scenarios  akin  to  the  one  described  in23  are  considered:  a 
stationary  A-elemcnt  transmitter  array  is  located  above  a  perfectly  flat,  lossy  ground;  the  scene  is 
sparsely  populated  with  scatterers,  among  which  are  the  targets  of  interest  (i.e.,  one  or  more 
nonlinearly  loaded  metallic  antenna  elements)  positioned  either  above  or  below  the  ground 
surface  (Fig.  1).  The  transmitter  array  is  composed  of  infinitesimal  dipoles,  emitting  either 
single-tone  or  multi-tone  time-harmonic  fields.  The  steady-state,  multistatic  scattered  signals — at 
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the  harmonic  frequencies  of  the  excitation  frequencies — from  the  scene  are  captured  at 
observation  points  co-located  with  the  transmitters. 


Fig.  1  Standoff  sensing  of  nonlinearly  loaded  antenna  structure  in  a  half-space  environment 
in  the  presence  of  clutter 


2.1  Linear  and  Nonlinear  Scattering  Characterization 


The  incident  fields  at  each  antenna  structure  are  calculated  with  the  dyadic  half-space  (one-layer- 
medium)  Green's  function.24  The  responses,  or  induced  currents,  on  the  structure  are  deduced  by 
exploiting  a  hybrid  method-of-moments  and  equivalent  circuit  approach;  the  antenna  input 
admittance  Yant  and  Norton  equivalent  short-circuit  current  Isc  are  first  characterized  at  the 

relevant  frequency  points  using  the  half-space  mixed-potential  integral  equation  solver,  and 

o 

then — in  the  circuit  domain — a  harmonic  balance  technique  is  employed  for  obtaining  the 
terminal  voltage  response  Vant  of  the  antenna  connected  to  the  nonlinear  load.  In  essence,  after 

dividing  the  loaded  antenna  circuit  into  linear  and  nonlinear  sections  at  the  antenna  terminal 
reference  plane  and  applying  the  standard  nodal  current  law,  the  following  nonlinear  matrix 
system  can  be  established  to  solve  for  the  terminal  voltage  response: 


■<g\ 


(w 


+  KntKnt  ~Isc+H 


v4 


V  -c 

r  ant  c 


=  0, 


(1) 


which  contingents  on  the  suppositions 
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P  is  the  number  of  frequencies  ( o/p  e  Q!s ,  or  /;f  e  F,s )  employed  in  the  hannonic  balance 
technique;  Q  is  the  number  of  excitation  tones  (o“eQ",o  r  /“  e  Fa ,  and  Q“  cz  Q/; , 

Fa  a  F/: );  and  /sr  are  defined  using  the  coefficients  of  Equations  2  and  3;  and  Tmf  is  the 

antenna  admittance  matrix  covering  the  P  frequencies  of  interest.  The  i-v  characteristic  of  the 
nonlinear  load  is  i(t)  =  g  ( v  ( /  )  j ;  transformations  between  the  time  and  frequency  domains  are 

obtained  with  the  matrices  Tf<_t  and  Tti_f .  Also  note  that  the  nonlinear  load  is  biased  with  DC 
voltage  Vb  and  resistor  Rh  (the  bias  inductor  is  ignored  in  the  AC  analysis);  thus,  in  compact 

form,  c  is  a  constant  vector  with  all  elements  equal  to  Vh  /  Rh  and  H  is  non-zero  only  for  the 
first  element:  Hn  =  1 .  The  solution  to  (1)  is  computed  with  the  Newton-Raphson  method,  after 
deriving  the  Jacobian  of  the  matrix  system  and  properly  defining  Tf<_t ,  Tti_f ,  and  the  time 

domain  sampling  rate  and  interval.  (Additional  definitions  for  variables  related  to  Equation  1  are 
in  Appendix  A.)  Subsequently,  having  deduced  Vant ,  the  induced  currents  on  the  antenna 

structure  at  each  frequency  are  found  with  the  aforementioned  method-of-moments  solver,  and 
the  scattered  fields  are  propagated  to  the  observation  points  using  the  same  method  as  for 
calculating  the  incident  excitation  fields.  In  computing  the  scattered  signal  at  the  illumination 
frequencies  ( a>pp  e  Q"  ),  the  short-circuited  antenna  scattering  pattern  is  added  to  the  antenna  re¬ 
radiation  component  to  obtain  the  total  scattering  response. 

Note  that  for  the  calculation  of  the  dyadic  half-space  Green's  function  within  the  method-of- 
moments  solver,  exact  Sommerfeld  integrals  are  used.  For  the  scattering  and  propagation  routine, 
however,  in  order  to  increase  the  computational  efficiency,  second-order-accurate  asymptotic 
approximations  to  these  Sommerfeld  integrals  are  employed  instead:-  ’  these  asymptotic  forms 
(Appendix  B)  are  derived  by  applying  the  method  of  steepest  descents  to  the  integral  expressions 
and  provide  the  necessary  accuracy  to  capture  wave  interactions  at  the  near-ground  region, 
including  Norton  surface  wave  effects  — which  can  be  of  importance  for  near-grazing  half¬ 
space  propagation  scenarios,  depending  on  the  physical  parameters  of  the  problem  and  the 
locations  of  the  transceiver  array  and  targets  with  respect  to  the  ground  interface. 

2.2  Subspace-Based  Image  Formation 

After  the  detennination  of  the  linear  and  nonlinear  scattered  signal,  a  subspace-based  imaging 
technique  similar  to  the  time -reversal  multiple  signal  classification  method"  that  has  been 
developed  for  linear  radar  sensing  is  adapted  for  target  localization.  The  mechanisms  and 
formulations  of  the  algorithm  implemented  herein  are  established  as  follows.  Considering  the 
case  of  Mpoint-like  targets  at  positions  rf  ( m  e  {1,2,...,M}  ),  and  denoting  the  current  on  the 

n-th  (  n  e  {1, 2, . ..,  N}  )  transmitter — at  frequency  coaq  — as 

■/„(-%  < ) = L  ( <  >(«„(<)  *  ■ +  (< )  r +  ? ,  (< )  d  <*(<•-■ K ) ,  (4) 
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the  incident  signal  is 
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es,  Gjrfx,^) 


4 


(6) 


in  which  G,n  ( r  f ,  r“ ,  of  j  is  the  dyadic  half-space  Green's  function.  Upon  noting  the  reciprocity 
relation  G7;  ( rf ,  #•“ ,  o/p  j  =  G^n  ( #•“ ,  rf ,  cop  )  ,  the  scattered  field  at  frequency  copp  at  all  /V 
receiver  positions  can  then  be  cast  as 

£,«)eC!™  =^(rV%fflpqrx(£,(n"),ffl;),  (7) 


with  £(•)  as  the  operator  transforming  excitation  fields  into  radiation  currents  at  rf  .  It  is  seen 

from  Equation  7  that  the  scattered  signal  is  a  subspace  spanned  by  the  columns  of 

Q[rp ,ra ,a>f)  (that  is,  the  Green's  function  vectors),  and  the  scatterer  positions  are  encoded 

within  this  representation.  Transmitting  the  3 N  dipole  element  components  sequentially,  a 
multistatic  scattered  field  matrix  K  [g>p')  e  CiNx3N  can  be  obtained  as 


*«H£K)  *;«)  -  *r«y 


(8) 


Accordingly,  the  signal  subspace  can  be  identified  with  a  singular  value  decomposition 
procedure  (that  is,  K  [o/p  )  =  U  V  [(Dp)  ),  and  the  imaging  functional  over  the 

interrogation  domain  V  at  oof  can  be  constructed  as 


f  IN 


z  z  “/Wr^*(r">  r’(°p) 

k= 1  V/=Z+1 


2  A 


(9) 


with  the  assumption  that  the  signal  subspace  is  spanned  by  the  first  L  singular  vectors 
u,[a>p) — which  correspond  to  the  first  L  significant  singular  values,  and 

Qk [ra ,  r,  cof )  =  Qk [r a ,  r,  (opp  )  /  ||  Qk  ( ra ,  r,  oof  )  || ,  where  the  vector  Gk[ra  ,r,  oof)  — generated  as 
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the  asymptotic  expressions  from”  (Appendix  B)  for  computational  efficiency  purposes — is  the 
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k- th  column  of 


.  The  objective  function 


GSn(r,K,<o*)  ew(r,r“,^)  ... 

in  Equation  9  can  also  be  computed  and  summed  over  multiple  hannonics  to  enhance  the 
stability  of  the  algorithm. 

As  relevant  to  the  linear  scattering  response,  in  the  vector  field  case,  note  that  each  metallic  point 
target  would  give  rise  to  a  maximum  of  three  significant  singular  values  (and  associated  singular 
vectors),  and  if  the  point  targets  are  well-separated,  each  of  these  sets  of  singular  values/vectors 
would  have  a  one-to-one  correspondence  with  the  targets.  An  extended  target,  in  contrast,  could 
potentially  give  rise  to  more  than  three  significant  singular  values/vectors.  The  imaging 
functional  of  the  form  in  Equation  9  is  expected  to  work  well  for  both  point  and  extended  targets, 
since  an  extended  target  could  be  approximated  as  a  collection  of  point  radiators,  each  with  a 
field  pattern  defined  by  (■) .  To  identify  the  null  subspace  needed  in  the  imaging  algorithm, 

the  condition  3 N  >  3 M  — where  M  is  the  total  effective  number  of  individual  scatterers  in  the 
scene — must  be  met.  Additionally,  as  mentioned  in  other  works,34'35  note  that  care  must  be  taken 
in  choosing  the  parameter  L,  especially  in  the  presence  of  noise. 

For  the  nonlinear  scattering  response,  it  is  seen  that  each  antenna  scatterer — regardless  of 
whether  it  can  be  treated  as  a  point  or  extended  target — has  only  one  significant  singular 
value/vector.  This  can  be  understood  by  noting  that,  in  essence,  the  scattered  field  at  each 
hannonic  frequency  co Pp  £  Of  is  found  by  exciting  the  antenna  tenninal  with  voltage  source 

Vant  [a>p  e  Q/;  j  — therefore,  the  receiver  array  is  capturing  the  antenna  radiation  pattern  and  not 

the  direct  scattering  response.  As  a  consequence,  without  a  priori  knowledge  of  the  hannonic 
radiation  pattern,  results  from  the  imaging  functional  in  Equation  9  may  be  degraded  in  some 
cases  for  the  extended  target  scenario.  To  improve  the  imaging  quality,  the  following  composite 
imaging  functional  can  be  used: 

where  the  linear  (or  ordinary)  scattering  imaging  result  at  excitation  frequency  of  is  combined 

with  the  nonlinear  (or  harmonic)  imaging  result  to  isolate  the  nonlinear  targets.  Assuming  there 
are  M  nonlinear  targets  in  the  scene,  note  that  for  the  construction  of  the  null  subspace  in  the 
hannonic  response  case,  the  condition  3 N  >  M  needs  to  be  satisfied. 


3.  Numerical  Experiment  Results 


The  framework  put  forth  in  Section  2  is  first  applied  for  analyzing  the  scattering  and  imaging 
responses  of  on-surface  and  subsurface  nonlinearly  loaded  antenna  targets  in  a  half-space 
environment;  then  the  complex-room  detection  problem  is  considered. 
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3.1  Sensing  in  Half-Space  Environment 


The  employed  sensing  geometry  is  similar  to  that  of  the  forward-looking  radar  system  described 
in:  ’  the  linear  transceiver  array  is  composed  of  N=  16  equally  spaced  elements  distributed 
over  a  2-m-wide  aperture  located  at  2  m  height.  In  the  forward-looking  direction,  this  sensing 
configuration  offers  resolution  primarily  in  the  down-range  and  cross-range  extents;  very  limited 
depth  resolution  in  the  vertical  extent  is  supported,  even  if  field  integration  over  multiple  moving 
apertures  is  utilized.  As  such,  the  only  localization  parameters  of  interest  are  the  down-range  and 
cross-range  positions  of  a  target.  In  the  first  example,  scenes  consisting  of  nonlinearly  loaded 
vertical  dipole  elements  (50  cm  in  length)  situated  on  the  ground  surface  are  considered.  The 
ground  soil  is  characterized  by  relative  dielectric  constant  er  =  4  and  conductivity  a  =  15mS/m. 

Figure  2a  shows  the  scattered  signal  spectrum  for  a  three-target  scenario  (M=  3)  in  which  each 
antenna  element  is  center-loaded  with  a  diode:  /(t)  =  Io  (exp(v(t)/v,)-l) ,  Io  =  lOnA, 

v,  =  26  mV ;  the  bias  network  has  parameters  =  1  kQ ,  V/;  =  1 V.  A  single-tone  excitation  (with 
transmitted  power  of  4  W)  at  f“  =  300  MHz  is  assumed,  resulting  in  a  scattered  field  spectrum 
composed  of  signals  at  multiples  of  the  excitation  frequency  af“  (a  e  {1,2,...}  ).  Note  that  only 
10  harmonic  frequency  components  ( |Q/J  |  =  P  =  10  )  are  computed  in  the  harmonic  balance 

solver.  In  this  example  (and  in  all  that  follow),  the  transmitter  array  is  centered  at  r  =  (-10  m,  2 
m,  2  m)  and  parallel  to  the  v-axis.  The  singular  value  spectrum  derived  from  the  multistatic 
response  at  600  MHz  is  illustrated  in  Fig.  3a,  which  shows  the  three  dominant  singular  values 
corresponding  to  the  targets.  Images  for  the  scene  constructed  with  Equation  9  using  the 
responses  at  frequencies  600  MHz,  900  MHz,  1500  MHz,  and  2100  MHz  are  displayed  in  Fig.  4, 
where  the  true  positions  of  the  targets  are  also  indicated.  Simulation  results  for  various  levels  of 
signal-to-noise  ratio  (SNR)  from  oo  (no  noise)  to  10  dB  are  obtained  and  examined;  the 
following  definition  for  SNR  is  observed: 


sNR=20iogl0(||jq®,qll/IU«) 


do 


in  which  N(copp  j  is  a  zero-mean  Gaussian  complex  noise  matrix,  and  ||  ■  ||F  denotes  the 
Frobenius  norm  operation.  The  interrogation  domain  here  is 

V  =  jr  =  (x,y,z)  e  R3 :  -5  m  <x,y<  5  m,  -1  m  <z  <  2mj,  and  the  xy  image  slice  that 
contains  the  maximum  response 


vp 

nl,max 


f 


=  max 

reO 


Z  X¥«i{r’(oP 

Kp= 2, 3, 5, 7 


(12) 


is  shown.  As  anticipated,  the  vertical  resolution  (in  z)  is  very  poor  in  this  case;  of  primary 
interest  is  the  imaging  performance  in  the  range  and  angular  directions.  (Throughout  this  work, 
the  parameter  L  in  the  imaging  functional  is  chosen  by  first  plotting  the  singular  value  spectrum 
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in  logarithmic  scale  and  then  identifying  the  “knee”  point  in  the  graph,  which  is  taken  as  the 
demarcation  between  the  signal  and  null  subspaces.) 


(a) 


300  600 

Frequency  [MHz| 

(b) 


900 


Fig.  2  Monostatic  scattered  electric  field  strength  for  center  array  element:  a)  Single-tone 
excitation  at  300  MHz;  three  diode-loaded  targets  and  b)  Two-tone  excitation  at 
300  MHz  and  305  MHz;  two  varactor-loaded  targets 
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Fig.  3  Normalized  singular  value  spectrum:  a)  Three  diode-loaded  targets  and  b)  Two 
varactor-loaded  targets 


9 


Fig.  4  Three  diode-loaded  vertical  dipole  antennas;  single-tone  excitation  (  Q“  |  =  Q  =  1 )  at  300  MFlz 

(transmitted  power:  4  W).  |  =  P  =  10  .  Imaging  formation  frequencies:  600  MHz,  900  MHz, 

1500  MHz,  2100  MHz.  a)  SNR  =  oo  (no  noise);  b)  SNR  =  30  dB;  c)  SNR  =  20  dB;  and  d)  SNR  =  10  dB. 
Circled  numbers  denote  nonlinear  targets.  The  image  intensity  is  displayed  in  dB  scale. 

As  expected,  the  presence  of  the  noise  tends  to  degrade  the  resolution  in  both  the  range  and 
angular  directions;  the  range  resolution,  in  particular,  appears  to  be  much  more  drastically 
affected.  It  is  noted  that  under  noiseless  conditions,  near  perfect  localization  is  seen — that  is,  as 
applicable  to  this  example,  the  targets  appear  as  near  point  responses  in  the  image  for  SNR  =  oo 
(Fig.  4a).  Specifically,  the  results  not  only  show  that  the  three  targets  are  accurately  localized  but 
also  numerically  confirm  the  imaging  method's  potential  in  achieving  super-resolution — which 
has  been  identified  as  one  of  the  advantages  of  using  a  subspace-based  method  ’  . 


10 


It  should  be  emphasized  that  the  reason  for  adopting  the  particular  aspect-limited  aperture  with 
the  prescribed  array  parameters  is  to  conform  to  the  operational  constraints  of  a  ground-based, 
forward-looking  system.  The  application  of  a  larger  linear  aperture  (or  a  full  aperture)  may 
improve  the  imaging  performance;  for  example,  as  applied  to  the  traditional  linear  scattering 
response,  full  circular  arrays  have  been  considered  in,  where  only  2D  problems  are  treated. 
From  the  implementation  perspective,  however,  large  aperture-  and  circular  aperture-based 
sensing  configurations — commonly  used  in  other  studies — are  not  practical  for  the  current 
application. 

To  investigate  the  decay  rates  of  the  harmonic  signal  intensities  in  terms  of  sensing  distance,  the 
scattered  fields  are  computed  and  plotted  as  a  function  of  standoff  range,  p  ,  or  the  radial 

separation  between  the  observation  point  and  a  target.  As  a  generalization,  with  small-signal 
analysis,  it  is  demonstrated  that — for  the  backscattering-based  sensing  scenario  in  which  the 
interaction  points  are  positioned  close  to  the  ground  surface — the  amplitude  of  each  harmonic 
order  of  the  scattered  field,  in  the  asymptotic  propagation  regime,  approximately  depends  on  the 
distance  as  oc  p  2<a"'1 ;  this  behavior  is  consistent  with  theoretical  results  derived  using  analytical 
models  for  the  nonlinear  circuit  response  and  field  propagation."  Thus,  the  nonlinear  scattered 
components  experience  faster  decay  rates  as  compared  to  the  conventional  oc  p  4  factor 
encountered  in  linear  radar  analysis;  this  propagation  characteristic  has  to  be  considered  in 
assessing  the  link  budget  in  the  design  of  the  sensing  system. 

The  scattering  and  imaging  results  for  a  2-target  scenario  are  displayed  in  Figs.  2b  and  5.  In  this 
case,  each  vertical  dipole  antenna  is  center-loaded  with  a  varactor: 

i(t)  =  Cj0  (l  +  v(t)/vw)  ■dv{t)fdl ,  Cj0  =  2pF,  vw=0.7V,  y  =  0.5  ;  the  bias  network  has 
parameters  Rh  =  1  kQ, ,  Vb  =  0. 1 V.  A  2-tone  excitation  (with  transmitted  power  of  1 .7  W)  at 
f“  =  300  MHz  and  =  305  MHz  is  considered,  which  results  in  a  scattering  spectrum 
consisting  of  responses  at  af“  +  bf“  ,  where  a  and  b  are  integers.  Note  that  only  12  harmonic 
frequency  components  ( |Q^  |  =  P  =  12 ,  |a|  + 1/>|  <  3  )  are  computed  in  the  harmonic  balance  solver. 

The  singular  value  spectrum  at  600  MHz  is  shown  in  Fig.  3b.  The  image  is  generated  using  the 
multistatic  scattered  signal  at  frequencies  600  MHz,  605  MHz,  and  610  MHz.  The  decay  rate  of 
the  harmonic  scattered  components  in  this  case  is  seen  to  be  of  the  form  oc  p1( H+l/’l+1) . 
Consequently,  the  relative  strengths  of  the  various  spectral  components  compared  to  each  other 
are  distance  dependent:  for  example,  although  the  5  MHz  (a  =  -  l,  b=  1)  component  is  weaker 
than  the  295  MHz  (a  =  2,  b  =  -1)  and  310  MHz  (a  =  —1,  b  =  2)  components  in  the  near  range 
(p  «  10  m),  it  may  become  the  stronger  component  at  longer  ranges.  The  main  gist  of  these 

observations  can  be  generalized  to  the  multi-tone  excitation  configuration. 
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Fig.  5  Two  varactor-loaded  vertical  dipole  antennas;  double-tone  excitation  (  Q“  |  =  Q  =  2  )  at  300  MHz  and 

305  MHz  (transmitted  power:  1.7  W).  |  =  P  =  12  .  Imaging  formation  frequencies:  600  MHz, 

605  MHz,  610  MHz.  a)  SNR  =  oo  (no  noise);  b)  SNR  =  30  dB;  c)  SNR  =  20  dB;  and  d)  SNR  =  10  dB. 


The  scenarios  considered  in  Fig.  6  are  similar  to  those  in  Figs.  4  and  5,  with  the  only  difference 
being  that  one  of  the  antenna  elements  is  treated  as  a  purely  linear  clutter  object;  in  Fig.  6a, 
scatterers  1  and  3  are  diode-loaded,  whereas  scatterer  2  has  no  nonlinear  loading  (i.e.,  in  this 
case,  the  antenna  is  short-circuited  at  its  feeding  terminal);  in  Fig.  6b,  scatterer  1  is  varactor- 
loaded  whereas  scatterer  2  is  short-circuited.  The  images  show  that  the  nonlinear  targets  are  all 
identified.  (For  space  considerations,  henceforth  only  the  results  for  SNR  =  20  dB  are  included 
here.)  It  should  be  mentioned  that,  for  the  simulation  parameters  and  the  set  of  inter-target 
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separation  distances  considered  in  this  example,  inter-scatterer  coupling  effects  are  negligible; 
for  smaller  separation  distances,  however,  it  is  possible  that  these  coupling  effects  could  become 
important — for  instance,  the  re-radiation  from  the  nonlinear  scatterer  could  induce  harmonic 
currents  on  the  linear  clutter  object,  which  in  turn  would  radiate  responses  at  the  harmonic 
frequencies  even  though  it  inherently  behaves  as  a  purely  linear  scatterer. 


Fig.  6  Same  scenarios  as  those  in  Figs.  4  and  5,  except  scatterer  2  is  a  linear  clutter  object.  SNR  =  20  dB 


In  the  next  example,  scattering  from  horizontal  dipole  antennas  (50  cm  in  length)  situated 
directly  on  top  of  the  ground  surface  is  simulated.  For  the  scenes  corresponding  to  the  images  in 
Fig.  7,  all  the  targets  are  nonlinear  except  for  scatterer  2  in  Fig.  7b,  which  is  short-circuited. 
Nonlinearly  loaded  targets  are  center-loaded  with  a  diode  that  has  the  same  parameters  as  those 
in  the  previous  cases.  It  is  seen  that  the  radiation  modes  from  the  horizontal  dipole  antenna  are 
not  as  well  approximated  by  G,q  (•)  as  those  from  the  previously  considered  vertical  dipoles; 

thus,  the  imaging  ansatz  in  Equation  9  is  applied  to  improve  the  localization  performance. 
Specifically,  single-tone  excitation  at  300  MHz  is  assumed,  and  the  fundamental  and  second 
harmonic  returns  are  captured  at  the  receiver  array.  Once  again,  the  “spreading”  of  the  target 
image  response  in  the  range  direction  is  primarily  due  to  the  presence  of  noise. 
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Fig.  7  a)  Two  diode-loaded  horizontal  dipole  antennas  and  b)  One  diode-loaded  horizontal  dipole  antenna 

(scatterer  1)  and  one  linear  clutter  object  (scatterer  2).  Single-tone  excitation  ( |  =  Q  =  1 )  at  300  MHz 

(transmitted  power:  4  W).  |Q^  |  =  P  =  10  .  Imaging  formation  frequencies:  300  MHz,  600  MHz,  and  SNR 
=  20  dB. 


In  the  last  example,  the  nonlinear  targets  are  embedded  in  a  clutter  field:  the  2  diode-loaded 
horizontal  dipole  antenna  targets  (scatterers  1  and  4)  are  either  located  on  the  ground  surface 
(Fig.  8a)  or  buried  in  the  ground  at  3  cm  depth  (Fig.  8b);  the  linear  surface  clutter  objects  consist 
of  small  metallic  cylinders  and  short-circuited  dipoles  of  different  orientations.  The  imaging 
functional  in  equation  9  is  applied  with  single -tone  excitation  and  fundamental  scattering 
observation  at  300  MHz,  and  harmonic  scattering  observation  at  1200  MHz.  The  images  indicate 
that  the  nonlinear  targets  are  detected  with  good  angular  resolution — albeit  the  range  resolution 
is  degraded  because  of  the  addition  of  noise  to  the  scattering  response. 
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Fig.  8  Two  diode-loaded  horizontal  dipole  antennas  (scatterers  1  and  4)  in  a  clutter  field:  a)  On-surface  and 

b)  Buried.  Scatterer  2  is  a  vertical  dipole;  scatterer  3  is  a  horizontal  dipole;  scatterers  5  and  6  are  metallic 

cylinders.  Single -tone  excitation  (|Q“|  =  Q  =  1)  at  300  MHz  (transmitted  power:  4  W).  Q/7  =  P  =  10  . 

Imaging  formation  frequencies:  300  MHz,  1200  MHz,  and  SNR  =  20  dB. 


3.2  Sensing  in  Complex-Room  Environment 

Localization  of  a  nonlinearly  loaded  antenna  target  is  considered  with  the  sensing  configuration 
shown  in  Figs.  9a,  10a,  and  1  la.  The  electrical  properties  of  all  the  structures  and  objects  (walls, 
ceiling,  floor,  furniture,  humans,  etc.)  found  within  the  complex  room  scene  are  described  in;40 
the  transceiver  array  is  similar  to  the  one  used  in  Section  3.1,  except  that  here  it  is  centered  at  r  = 
(6.5  m,  0  m,  2  m);  and  the  target — centered  at  1  m  height — is  the  same  diode-loaded  vertical 
dipole  antenna  that  was  previously  noted.  The  procedure  for  scattering  characterization  closely 
follows  the  one  outlined  in  Section  2.1,  with  field  propagation  to  and  from  the  target  computed 
using  a  finite-difference  time-domain  solver.40  As  the  exact  Green’s  function  for  this 
environment  is  complicated  and  cannot  be  readily  derived,  one  key  approximation  is  made  to 
facilitate  the  analysis;  for  both  the  integral  equation  solver  and  the  image  generation,  the  half¬ 
space  Green's  function  is  assumed.  The  imaging  results  for  three  target  positions  are  displayed  in 
Figs.  9b,  10b,  and  1  lb:  in  all  cases,  a  single-tone  excitation  at  f“  =  300  MHz  is  employed,  and 

the  images  are  constructed  using  the  responses  at  600  MHz.  It  is  seen  that  while  the  cross-range 
resolution  is  satisfactory,  the  range  resolution  is  rather  poor.  The  degradation  in  range 
resolution — as  well  as  in  range  positioning  accuracy — is  primarily  due  to  the  presence  of  the 
front  wall  and  the  in-room  clutter  in  the  vicinity  of  the  target — the  effects  of  both  of  which  are 
not  taken  into  account  in  the  imaging  construction,  given  that  only  the  half-space  Green's 
function  is  applied  within  the  imaging  functional  Equation  9. 
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Target 


x  [m] 


(b) 

Fig.  9  a)  Diode-loaded  dipole  antenna  in  a  complex-room  environment — 

position  1  and  b)  Imaging  result  with  single-tone  excitation  ( |  =  Q  =  1 )  at 

300  MHz  (transmitted  power:  4  W).  |  =  P  =  10  .  Imaging  formation 

frequency:  600  MHz,  SNR  =  20  dB. 
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Fig.  10  a)  Diode-loaded  dipole  antenna  in  a  complex-room  environment — 

position  2  and  b)  Imaging  result  with  single-tone  excitation  ( |  =  Q  =  1 )  at 

300  MHz  (transmitted  power:  4  W).  |  =  P  =  10  .  Imaging  formation 

frequency:  600  MHz,  SNR  =  20  dB. 
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Target 


(a) 


Fig.  11  a)  Diode-loaded  dipole  antenna  in  a  complex-room  environment — 

position  3  and  b)  Imaging  result  with  single-tone  excitation  ( |  =  Q  =  1 )  at 

300  MHz  (transmitted  power:  4  W).  |  =  P  =  10  .  Imaging  formation 

frequency:  600  MHz,  SNR  =  20  dB. 
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4.  Conclusions 


To  facilitate  the  development  of  a  complementary  nonlinear  detection  capability  for  an  existing 
near-ground,  forward-looking  linear  imaging  radar  system,  a  numerical  solver  is  presented  in  this 
work  for  characterizing  the  multiview,  multistatic  RF  harmonic  response  of  nonlinearly  loaded 
antenna  targets  under  tone  excitation — viz.,  a  half-space  method-of-moments  approach  is  applied 
for  deriving  the  equivalent  circuit  of  the  antenna  structure;  the  harmonic  balance  technique  is 
employed  in  solving  for  the  steady-state  response  of  the  nonlinear  system;  and  asymptotic 
methods  are  adapted  when  appropriate  in  the  propagation  and  scattering  routines  to  characterize 
the  re-radiation  modes  of  the  scatterer.  Additionally,  a  subspace-based  direct  imaging  technique 
is  developed  and  demonstrated  for  near-field  location  estimation  of  on-surface  and  buried 
targets.  An  extensive  set  of  numerical  experiments  has  been  carried  out,  and  scattering  and 
imaging  results  have  been  derived  for  sensing  scenarios  with  and  without  linear  clutter.  The 
imaging  responses  indicate  that  by  using  an  aspect-limited  transceiver  aperture  consistent  with 
that  of  the  aforementioned  forward-looking  radar,  accurate  target  localization  can  be  obtained.  In 
sum,  the  harmonic  detection  technique  is  observed  to  be  a  viable  method  for  isolating 
nonlinearly  loaded  targets  from  surrounding  linear  surface  clutter  in  a  sparse  scene;  however,  the 
imaging  performance  is  determined  to  be  rather  noise  sensitive — with  the  down-range  resolution 
being  much  more  vulnerable  to  noise  effects  than  the  angular  resolution;  the  practicality  of 
nonlinear  sensing  is  also  limited  to  the  near  range  due  to  the  higher-order  attenuation  rates  of  the 
harmonic  backscattered  fields.  Although  the  focus  of  the  study  is  on  the  half-space  problem, 
target  localization  and  imaging  in  a  complex-room  environment  are  also  considered.  It  is  seen 
that  the  use  of  an  approximate  Green's  function  (e.g.,  the  half-space  formulation)  in  lieu  of  an 
exact  formulation  in  the  imaging  objective  function  can — as  expected — lead  to  degradation  in 
the  range  resolution  and  positioning  accuracy,  whereas  the  angular  resolution  and  positioning 
accuracy  are  less  affected. 

Note  that  the  electromagnetic  scattering  and  imaging  approaches  presented  herein  can  be 
extended  for  the  treatment  of  wideband  sensing  scenarios.  For  example,  the  targets  can  be 
excited  by  multiple,  separate  sets  of  tones  in  a  sequential  manner,  and  then  the  imaging 
functional  is  integrated  over  all  harmonic  response  sets  to  enhance  the  detection  performance.  A 
moving  transceiver  aperture  configuration  could  also  provide  multiple  independent  scattering 
responses  needed  to  further  improve  the  accuracy  of  the  localization  scheme.  These  issues  will 
be  investigated  in  future  studies. 
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The  following  definitions  are  needed  for  Equation  1 : 

V  e  M(2/>+1)xl  =  \v  V  V  V  V  T-  tA-11 

'  ant  c  ^  L  ant, DC  r  ant, l  r  ant, 2  •••  r  ant,2P-\  r  ant, IP  J  ’  V*  1  / 

the  elements  of  /  e  M(2/mh  are  only  nonzero  at  the  excitation  frequency  components®"  : 


Jsc,2q-\  I sc,2q  _  ’ 


(A-2) 


the  antenna  admittance  matrix  is 


Y  e  ]g>(2-p+1M2-p+1)  = 

ant 
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0  G(0)f)  B(cof) 

0  -B(o)f)  G(cof) 


G(®f)  5(®f) 

-5(®f)  G(®f) 


G{copP)  B(a>p) 
-B(a>pp)  G(a>p) 


with  G(®p  =  «(K„,„(fflp)  and  = 

Denoting  the  temporal  sampling  rate  as  St  and  the  total  number  of  time  steps  as  Nt ,  the 
transformation  matrices  can  be  shown  to  be 


TV/  e  . 


,JV,x(2P+l) 


(A-3) 
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and  Tf<_t  e  R(2P+1)xAr'  =  ( 


tt  t  \ 1  Tt 

±t^f  ±t^-f  /  At<-f  • 


(A-4) 


In  accordance  with  the  i-v  characteristics  employed  in  Section  3,  the  Jacobian  of  the  left  hand 
side  of  Equation  1  with  respect  to  Vant  is  derived  as 


•7(r,,)  =  7}«dia  g 


di(t ) 
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dv(t) 

KmM 

\/t=ntSt  j 

t.  ,  +  r  +—h 


t4-f  ant 
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for  a  diode,  and 
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j(yM)=T} 
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for  a  varactor,  in  which  nt  e  {0,1,...,  JV,  -1}  ,  c(v(t))  =  Cj0  (l+  v(*)/vw)  ?  ,  and 
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Intentionally  Left  Blank. 


26 


27 


Closed-form  analytical  expressions  for  G,  (•)  are  derived  from  exact  Sommerfeld  integrals 

using  the  method  of  steepest  descents.  Based  on  the  near-ground  sensing  geometry  presented 
herein,  using  the  notational  convention  from,26  G^  (•)  =  G(m  (•)  =  G f  (■)  +  G  (•)  for  on-surface 

targets  and  G,n  (•)  =  Gm  (•)  =  Gt  (•)  for  buried  targets,  where  Gfs  (•)  is  the  free-space  wave 

component.  By  following  the  derivation  outlined  in,26  for  an  infinitesimal  electric  current 
element  located  on  the  z-axis  at  z'  and  an  observation  point  at  r  =  (p,  <j>,  z ) ,  it  can  be  shown  that 

the  second-order-accurate  ground  wave  components  Gr  (•)  and  G,  (•)  are  given  by 
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where  the  saddle  point  is  defined  as  cos  ws  =  (z  +  z')  /  r  for  Gr  (•)  and  cos  ws  «  z  /  p  for  £,(•), 
and  the  following  ground  reflection  and  transmission  coefficient-related  expressions  are  needed: 

mr  ±  (w)  =  r  ( w)  ±  r+  ( w)  cos(2 <j>) ;  (B-4) 

mt  ±  (w)  =  t+(w)±t_{w)  cos(20) ;  (B-5) 
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